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We present a theory of time-dependent tunneling between a metal and a partially spin-polarized 
two-dimensional electron system (2DES) . We find that the leakage current which flows to screen an 
electric field between the metal and the 2DES is the sum of two exponential contributions whose 
relative weights depend on spin-dependent tunneling conductances, on quantum corrections to the 
electrostatic capacitance of the tunnel junction, and on the rate at which the 2DES spin-polarization 
approaches equilibrium. For high-mobility and homogeneous 2DES's at Landau level filling factor 
i^ = f , we predict a ratio of the fast and slow leakage rates equal to {2K + 1)^ where K is the number 
of reversed spins in the skyrmionic elementary charged excitations. 
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There has been renewed |l[] interest recently in non- 
equilibrium spin accumulation [0 due to electronic trans- 
port in spin-polarized electron systems, in part because 
these accumulations are important in giant magnetore- 
sistance p]. In this paper we address spin-accumulation 
in the tunneling current between a metal and a two- 
dimensional electron system (2DES) in the quantum Hall 
regime. Our work is motivated in part by recent exper- 
iments W which have discovered unexplained two-rate 
leakage currents for such tunnel junctions, and also in 
part by the long M spin-relaxation times of 2DES's, es- 
pecially long in the quantum Hall regime P-H|. We find 
that spin accumulation depends subtly on the interplay of 
spin-dependent tunneling conductances, thermodynamic 
densities-of-states, and spin relaxation rates. According 
to our theory, the double-rate current found in experi- 
ment signals sizable quantum corrections to the effective 
capacitance of the junction and spin-relaxation processes 
which depart from those in conventional metals, both 
expected in the quantum Hall regime. Gapped quan- 
tum Hall states lead to rapid variation of the chemical 
potential ^ with density, and the presence of Skyrmion 
elementary charged excitations |10| requires an unusual 
spin-relaxation process. We predict that, for homoge- 
neous 2DES and thin tunneling barriers, the ratio of fast 
and slow relaxation rates dX v = 1 will equal {2K + 1^ 
, where K is the number of reversed spins in the lowest 
energy Skyrmion quasiparticles. 

We start from the following phenomenological linear 
response equations which we believe to be broadly appli- 
cable for tunneling between a metal and a 2DES. 






(1) 



Here Q^ and /Xo- are the spin-cr particle-number and 
chemical potential in the 2DES. These equations express 
spin-partitioning of the tunneling current; the assump- 



tion of separate chemical potentials for the 2DES spin 
subsystems is valid when the spin-relaxation time is much 
longer than other characteristic scattering times in the 
2DES. In these equations we have placed the zero of en- 
ergy at the chemical potential of the metal. The two 
terms on the right hand side of Eqs. ||] account respec- 
tively for tunneling across the junction and relaxation of 
the 2DES spin subsystems toward mutual equilibrium. 
A closed description of electron transport in the system 
requires, in addition to the above conductance equations, 
which relate currents to chemical potential differences, a 
set of capacitance equations which relate these chemical 
potentials to accumulated charges: 

A.T = -^o + (C7-i)TTQT + (^"')nQi (2) 

Iii^-Vo + {C-')i^Q^ + {C-')liQi (3) 

Here Vq represents the electrostatic contribution from 
charges external to the 2DES. Elements of the inverse 
capacitance (C~^) matrix have an electrostatic contribu- 
tion proportional to the width of the tunnel barrier and 
a quantum 'chemical potential' contribution due to the 
Fermi statistics and correlations of electrons in the 2DES: 
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(4) 



where Cg = Ae/{A'Ke^d) is the electrostatic capacitance 
of the junction, A is the cross-sectional area of the two- 
dimensional electron system, e is the dielectric constant 
of the host semiconductor, d is the distance between the 
metallic electrode and the 2DES, and jla is the spin-a 
chemical potential of the 2DES relative to its electric 
subband energy. The notation above is motivated by the 
relationship between -Fg-.o-' and Fermi liquid theory inter- 
action parameters. In the commonly employed Hartree 
mean-field approximation. 
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(5) 



where D° is the non-interacting 2DES density-of-states 
per spin. In general dfi^/dr]^i is non-zero for a ^ a' 
because of electronic correlations. 

These equations can be used to describe various time- 
dependent or ac transport experiments; we apply them 
here to the situation studied recently Q by Chan et al. in 
which a chemical potential difference across the junction 
is created by external charges and the leakage current 
which reestablishes equilibrium is measured as a function 
of time. In an obvious matrix notation we rewrite the 
conductance and capacitance equations as Q = G/z and 
/i = — Vo -I- C~-'-Q. Eliminating the chemical potentials 
using the capacitor equations yields a set of two coupled 
first-order inhomogeneous linear differential equations for 
the time-dependent spin-up and spin-down charges in 
the 2DES. Solving these with the boundary condition 
Q^{t = {)) — Qi{t = Q) = yields for the spin-dependent 
currents into the 2DES: 

Qa{t) = Ia,+ cxp(-i/r+) + la,- exp(-i/r-) (6) 

where t^ and tZ , generalized RC relaxation rates , are 
the eigenvalues of A = GC^-*-. We have obtained the 
following explicit expressions for /|.±: 



(7) 
(8) 



The corresponding expressions for I^ ± are obtained by 
interchanging spin labels. In Fig. [Q we show results 
obtained for the dependence of leakage current on 2DES 
chemical potential near Landau level filling factor ly = 1 
which follow from a non-interacting Skyrmion model of 
the 2DES described below. These replicate all major 
features found in experiment M. The peak in both fast 
and slow relaxation rates is due to the sharp decrease in 
capacitance as the 1^ = 1 incompressible quantum Hall 
state is approached. The leakage current is dominated 
by the slow process, except in a narrow range very close 
to :/ = 1 where the fast process takes over. The origin of 
this crossover is explained below. 

Similar results can be obtained for the instantaneous 
chemical potentials of the spin-up and spin-down subsys- 
tems: 




fJ-a{t) = -Vo +^ Ha,s{'^ - ex-p{-t/Ts)) 



(9) 



where fia,± — X^o-' ^ct<^"^o-',±'''±- We note that ii^,^ + 
fj-a,- = Vo; current flows until the electrochemical po- 
tential change for each spin cancels the electric potential 
from external charges. The two spin subsystems are in 
equilibrium at both the beginning and the end of the re- 
laxation process, but are, in general, out of equilibrium 



at intermediate times. The non-equilibrium spin accu- 
mulation /i|(i) - ^x(t) = (^|^_ - /X|,_)(exp(-i/r_) - 
exp(— t/T_(.)). What is readily separated in experiment 
are the fast and slow relaxation contributions to the 
current, not the spin subsystem contributions. Never- 
theless, we see that non-equilibrium spin accumulations 
occur system between times r_ and r_|_ whenever both 
contributions are present. In Fig. l2l we plot, time- 
dependent chemical potentials calculated for the non- 
interacting Skrymion model at /x = 0.05(e^/£) . The 
chemical potentials change quickly on a time scale r"*" 
and fi-^ overshoots Vb; on a much longer time scale the 
two chemical potentials approach Vq from opposite sides. 
We note from Fig. ||l| that the slow leakage current pro- 
cess actually dominates the capacitance at this value of 
fi: the naive view that the fast process is current flow to 
the 2DES while the slow process is spin-equilibration is 
incorrect. 

Before turning to the quantum Hall regime, where non 
equilibrium spin-accumulations are large, it is instructive 
to examine several limits for which spin-accumulation 



does 



not occur. For 

= (Gt + Gx)/Cg, 



Fa.a'Ca « 1 
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-1 



0, /, 
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we find that 
= GcrVo and 



/cr.- -^ 0. In this limit, which usually holds for metallic 
electrodes, spin-independent electrostatic contributions 
dominate electrochemical potential changes; no non- 
equilibrium spin-accumulation occurs because the spin 
subsystems are not driven from equilibrium by electro- 
static potentials. For strong tunnel barriers (Gct << Gs), 
on the other hand, t^^ = Gs(F|,| + Fi^i - 2i^-f,x), 
tZ = (G| -I- Gx)(F|_|F;,| - F^2^j^)/(i^-f,| -I- Fi^i ~ SF^^^) 
and that the fast relaxation current /+ = /|,+ -|-/|_+ — 0. 
For this limit spin accumulation does not occur be- 
cause the relaxation processes are fast enough to main- 
tain instantaneous equilibrium. Unlike the electrostatic- 
dominance case discussed first, the leakage current flows 
at the slow rate tZ ■ A third more subtle limit in which 
spin-accumulation does not occur applies to Fermi gas 
2DES's in which we may ignore correlation contributions 
to the chemical potential and the commonly adopted 
forms Ga — cD^, Gg = CgD^Di hold. These expres- 
sions result from golden rule estimates of quasiparticle 
tunneling and spin-flip transition rates respectively , c is 
a constant which declines exponentially with the thick- 
ness of the tunneling barrier and Cg is a constant de- 
pendent on spin-orbit scattering strength in the 2DES. 
For this model we find that T_jr^ = c(I/Gg + D^ + Di), 
tZ^ = {c+Cs{D^+Di)/Cg and all the weight is in the fast 
leakage current. No spin-accumulation occurs because 
the spin subsystems are not coupled by interactions and 
the ratio of tunneling conductances equals the ratio of 
the rates at which the chemical potentials increase with 
density. Finally we mention the case in which the 2DES 
is paramagnetic to which we return below. For G| = G| 
and i^ii = Fidi symmetry forbids spin accumulations. 



An explicit calculation finds no weight for the slow leak- 
age current and the rate ratio 

r_ _ GT(2/gg + -FT,T+na) 
r+ (GT+2G,)(Ft,T-FT,x)- 

None of these limits apply throughout the quantum 
Hall regime. Near integer Landau level filling factors, 
Fermi statistics and correlations in the 2DES, not electro- 
statics, dominate the electrochemical potential changes 
with density H. Equilibrium electronic states contain 
prt complex Skyrmion quasiparticlcs whose formation 
from the fully spin-polarized v ~ I ground state can- 
not be achieved by a single-particle process. Spin- 
equilibration will therefore be slow [|lO[. The two spin 
systems are intricately coupled so that the Fermi gas limit 
does not apply. Furthermore, the 2DES will generally be 
strongly spin-polarized. 

A simple model of the 2DES valid at low temper- 
ature for ly near one, is obtained by ignoring interac- 
tions between Skyrmions. We obtain the following grand- 
canonical ensemble expressions for the occupation prob- 
abilities of the N^ = A/(2Tr£'^) Skyrmion quasielectron 
and quasihole states with K excess reversed spins: 



n-Kh = f{eK + {K + 1)^1 - Km). 



Here /(e) = (exp(e/fcBT) + 1)-^ is a Fermi factor 0, 
ek is the energy of a Skyrmion quasiparticle, (27r£^)~^ 
is the density of a full Landau level, and we have cho- 
sen the zero of energy so that quasielectron and quasi- 
hole skyrmion states have the same energy ||l^. When 
the spin subsystems are in equilibrium (/i| ^ ji^ = fi) 
we can use Eqs. (12) to calculate the chemical potential, 
given the Landau level filling factor. Eqs. (n2|) express 
the property that the K-th quasielectron Skyrmion is 
formed by adding K -\- 1 spin-down electrons and remov- 
ing K spin-up electrons from to the v = 1 ground state. 
For non-interacting electrons only A' = quasiparticlcs 
occur; for typical 2DES's, on the other hand, the lowest 
energy quasiparticlcs have AT = 3 [|l3|. From Eqs. ( p^ ) 
we calculate the following thermodynamic densities of 
states, I?a-,CT' = dn^/dfi^i for the coupled spin systems: 



and A|^| occur in the ratio (AT + 1)^ : K'^ : A'(A: -|- 1) 
= 16 : 9 : 12. This contrasts with the non-interacting 
electron case for which Af^f is much larger than F^^^ and 
(10) -^T,i vanishes. (The same ratios apply for z^ < 1 with 
inverted spin- indices.) For u = 1 the low-temperatures 
equilibrium state charge is added to the v = I state in 
the form of AT = 3 skyrmions, i.e. for each 4 up-spins 
added to the 2DES, three down spins are removed. It 
follows that the time-integrated spin-up and spin-down 
leakage currents are approximately equal in magnitude 
and opposite in sign. 

In contrast to the partitioning of total leakage charge 
between spins , which is determined purely by thermo- 
dynamic considerations, the observable partitioning be- 
tween fast and slow components is difhcult to under- 
stand intuitively in the general case. Simplification oc- 
curs, however, when v = 1. It follows from particle- hole 
symmetry [|2[, that the symmetries for leakage into a 
paramagnetic 2DES also hold at i^ = 1, explaining the 
vanishing weight in the slow leakage current channel seen 
in Fig. ||l| as the filling factor approaches one. The ra- 
tio of fast to slow leakage rates, (Eq. ( [l(]| ) depends only 
on thermodynamic quantities provided that the tunnel- 
ing barrier is thin enough that Go- >> Gs. Then, pro- 
vided that the temperature is sufhciently low that quan- 
tum terms dominate the inverse capacitance we find that 
t^/t+ = (2 AT + 1)^ = 49, in rough agreement with the 
findings of Chan et al. M. We ascribe discrepancies to 
the inhomogeneity present at integer filling factors in all 
current samples. 
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27re^D^^^ =J2[iK + ifAieK + ^i) + K^A{eK ~ ^i)] 

K 
K 

where A(a:) = sech^{x/2kBT)/'ikBT. In Fig. [| we 
plot quantum contributions to the model's inverse capaci- 
tance, obtained by inverting this density-of-states matrix. 
The large peaks in capacitance near v — 1 are responsi- 
ble for the peaks in both leakage current rates. In the 
low-temperature limit, only the lowest energy AT — 3 
skyrmion will contribute so that for i/ > 1, A^j A^^j, 
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FIG. 1. Leakage current between a metallic electrode and 
a 2DES as a function of chemical potential near f — 1. The 
chemical potential is in e^/£ units, and the zero of energy 
is chosen so that/i = at i^ = 1. The solid and dotted 
lines show leakage rates in fast and slow channels respec- 
tively, while the long and short dashed lines show the ca- 
pacitance contributions, /+r+, and 7_r_, from fast and short 
rate channels respectively. These curves were calculated using 
the non-interacting Skyrmion model explained in the text to 
evaluate the quantum inverse capacitance contributions and 
a separation d — 51 between the metal and the 2DES. Ga 
and Gs were assumed to have the form 1000/Fcr,cr + 10 and 
l/(-^T,T-^J.,i) + 0-01 i'^ arbitrary units; the two terms repre- 
senting uniform system golden rule and inhomogeneity con- 
tributions respectively. The leakage rates are in units of C^^ 
times the conductance units and the capacitances per unit 
area are in units of £~^ . This figure is for ksT = 0.025e^/£. 
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FIG. 2. ficr + Vo in units of Vb as a function of time for 
the model of Fig. [|l| at /i = 0.056^/^". The maximum time 
illustrated is, 10~* units so that t/r^ large but f/r_ is still 
small at the largest times. 
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FIG. 3. Quantum contributions to the inverse capacitance 
for a 2DES near u = 1. The sohd lines show the majority 
and minority spin results for Fa.cr for a non-interacting model 
which contains only K = quasiparticles. For the Skrymion 
model, -Ft, I, -F^j and F^^^ are shown by dotted, short-dashed, 
and long-dashed lines respectively. 



